Intermediate Haldane phase in spin-2 quantum chains with uniaxial anisotropy 
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We provide evidence of an intermediate Haldane phase in a spin-2 quantum chain. By combining 
effective field theory and numerical approaches, we show that the phase diagram of the proposed 
model includes SO(5) Haldane, intermediate Haldane, and large-D phases. We determine the char- 
acteristic properties of these phases, including edge states, string order parameters, and degeneracies 
in the entanglement spectrum. The symmetries responsible for the degeneracy patterns observed in 
the entanglement spectrum are also discussed. 
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Introduction. Characterization of quantum phases be- 
yond Landau's symmetry breaking paradigm is an im- 
portant open problem in physics. Recently, for one- 
dimensional (ID) gapped phases, a classification scheme 
based on matrix product states (MPS) has been put 
forward "\ based on the fact that ground states of ID 
gapped Hamiltonians can be efficiently approximated by 
MPS'"'. This scheme complements the conventional 
approach that classifies relevant perturbations of fixed 
points, and sheds new light on some phases that have 
been extensively studied, e.g., the Haldane phase in 
integer-spin chains'. 

These developments may be helpful in order to clar- 
ify a controversial problem: the possibility of an inter- 
mediate Haldane (IH) phase (also called intermediate- 1? 
phase) in quantum spin chains. This problem originated 
from the study of a spin-2 Heisenberg chain with uniaxial 
anisotropy H = J2-j Sj ■ Sj+i+ DJ2j{S^)^, where D > 

(here S and S^ are the usual spin operators). For D — 0, 
the ground state is in the so-called Haldane phase. For 
i? — >■ oo, the ground state is in a large- Z? phase, close 
to a trivial product state with Sj=0 for all j. Regard- 
ing the phase diagram of this model, a field theoretical 
approach'" suggests a single phase transition between the 
Haldane and the large-D phases. On the contrary, for 
intermediate D, Oshikawa ' predicted that an IH phase 
may emerge between the two phases. He justified this by 
noticing that Sj = ±2 states are substantially suppressed 
by the D term but Sj — ±1 states are less affected, which 
may lead to the formation of an effective spin-1-like Hal- 
dane phase with residual S^ = ±1,0 states. Neverthe- 
less, whether such an IH phase really exists or not in 
the above spin-2 model (and its generalizations) remains 
unclear'""' ^ 

In this work, we study a spin-2 quantum chain for 
which we provide sharp evidence of the existence of an IH 
phase. As far as we know, our results provide the clearest 
evidence so far in favor of the existence of such a phase 
in quantum spin chains. Our search for the IH phase is 
guided by an effective field theory, which yields a qualita- 
tive phase diagram with three quantum phases, namely: 
S0(5) Haldane, IH, and large-I? phases. The field the- 
ory also determines characteristic features of the phases, 



e.g., edge states in open chains and string order parame- 
ters (SOPs) . We also determine numerically the phase 
diagram of the spin model, and find full agreement with 
the field-theory predictions. Moreover, we study the en- 
tanglement spectrum (ES)"' of the different phases, and 
confirm that the IH phase has a double degeneracy in the 
ES that is protected by symmetries, which distinguishes 
itself from the SO (5) Haldane phase with quadruple de- 
generacy and the large-D phase without protected degen- 
eracy. The spatial inversion, time reversal, and {Z2 x ^2)^ 
symmetries responsible for the robust degeneracy of the 
ES are also investigated. 

Model Hamiltonian and symmetries. In this work, we 
consider the spin-2 quantum chain 



H = Y.T.MSrs,+ir + Dj2{s[^ 



(1) 



for periodic boundary condition and in the thermody- 
namic limit. In our case, we take Ji = —^■,J2 = 

~Tio^ "^3 ~ §5' "^4 ~ el)' ^^'^ consider the region D > 0. 
As we shall see, this generalization of the usual spin-2 
Heisenberg chain has a number of important properties. 
First, let us identify the symmetries of Eq. (1), which 
will turn out to be very useful for our purposes. For 
D = 0, the Hamiltonian in Eq. (1) can be rewrit- 
ten as Hd=o - '2J2j[P2U,J + 1) + Pi{J,J + 1)], where 
PsriJiJ + 1) projects onto total spin-^T states of neigh- 
boring sites j and j + 1. This model has S0(5) symmetry 
and an MPS as its exact ground state'^'"'^'^. To identify 
the SO (5) symmetry, we work in the standard S^ basis 
\m) ijn = ±2, ±1,0) and define S0(5) Cartan generators 
ii2 ^ |2)(2| - I - 2)(-2| and L^^ = |1)(1| - | - 1)(-1|. 
By defining i^^ = ^(|2)(0| + |0)(-2| + h.c.) and i^s = 



V2 



(|1)(0| + |0)(-1| +h.c.), the S0(5) commutation rela- 
tions [L'^^ L'^''] = i {SacL'''' + SmL"-"^ - SadL''^ - SbcL"-"-) 
fix the ten generators L°-^ (1 < a < 6 < 5). For D = 0, 
the Hamiltonian commutes with all ten operators ^ L'f' 
and therefore has SO (5) symmetry. 

For I? > 0, this S0(5) symmetry is explicitly broken 
down to U(l)xU(l). In order to see this, we rewrite 
the uniaxial anisotropy as (S"^)^ = 4(L^^)^ + {L'^^Y. 
Thus, ^ LY and ^ L^'^ commute not only with each 



other but also with H, and therefore the model has 
U(l)xU(l) symmetry. Additionally, the Hamiltonian 
in Eq.(l) also has discrete symmetries, including spa- 
tial inversion, time reversal, and a set of Z2 symme- 
tries. These Z2 symmetries are a consequence of the 
invariance under global Z2 rotations [/"'' — e"^ for all 
L""^ . The Z2 operators form a {Z2 x Z2Y group, whose 
elements, without loss of generality, can be chosen as 
{1, C/12} X {1, [/15} X {1, ^734} X {1, [/35}. Here we re- 
mind that these Z2 operators preserve their form under 
the nonlocal Kennedy- Tasaki transformation and also 
generate a hidden {Z2 x ^2)^ symmetry in dual space^' . 
Field-theory treatment. Even though the model in Eq. 
(1) looks quite complicated, its effective field theory at 
low energy is very simple. As we shall explain soon, this 
is given by the following Hamiltonian density of five Ma- 
jorana fermions ^^ (a = 1, . . . , 5): 
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FIG. 1: (Color online) Qualitative phase diagram of the spin-2 
model in Eq. (1). The diagrams in the upper side correspond 
to the MPS structure describing the renormalization group 
fixed point for each phase: for low-D, ancillary spin-3/2 sin- 
glets are projected onto spin-2 subspaces at each physical site, 
thus providing an MPS for spin-2 physical particles in terms 
of matrices of dimension 4; for intermediate-D, ancillary spin- 
1/2 singlets, together with on-site spin-1/2 triplets S^ = 0, 
are projected onto spin-2 at each site; for large-D, no singlets 
are projected and the state is a product state. 



phase emerges, whose characteristics will be discussed be- 



U^e = -iv ^{i%d^£,% - S.ld.^S.1) - imi V Cr^I lo^- For D > D^^, the system enters the large-L" phase. 



4 
a— 3 



(2) 



where v and m^ are velocity and masses of the Majo- 
rana fermions, and marginal four-fermion interactions 
have been neglected. 

The strategy to derive Eq. (2) is to start from the 
S0(5) point Z? = 0, whose effective field theory is known 
to be of the form given by Eq.(2) with m\ = TO2 = ra^ < 
0^". In the continuum limit, the effect of the D term 
can be taken into account by using bosonization tech- 
niques. Following Ref.-"'-', we find that J2A-^]'^)'^ "^ 
wIdxj:L,eReL and E,(^f )' - ^g ^ dxYl^^^dl, 
where g < 0. By using [Slf = A{Lff + {iff, we ar- 
rive then at the expression given in Eq. (2). For D — >■ 0, 
we have mi — m^ ~ 4(7n2 — m^) ex D, but this relation 
does not hold for larger D due to renormalization ef- 
fects. Thus, the Majorana masses and the velocity v are 
treated as phenomenological parameters. However, the 
fact that only three independent masses appear in Eq. 
(2), which ensures 0(2) x 0(2) symmetry, is imposed by 
the U(l)xU(l) symmetry of Eq. (1) since 0(2)~U(1). 
Moreover, the {Z2 x ^2)^ symmetry of Eq. (1) is also 
revealed by the invariance of Eq. (2) under Z2 transfor- 
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Armed with this effective field theory description, we 
are now in position to sketch a phase diagram for the 
quantum spin chain in Eq. (1). When increasing D from 
to 00, we expect that the Majorana mass 7713 is al- 
ways negative in Eq. (2), while mi and 7712 change from 
negative to positive successively at two quantum critical 
points Dc^ and D^^ {Dci < Dc^). Both critical theo- 
ries at Dci and I?C2 have two massless Majoranas, and 
thus are equivalent to conformal field theories with cen- 
tral charge c = | x 2 = 1. For < D < Dc^, we call the 
phase 'S0(5) Haldane phase', since its physics is captured 
by the S0(5) point D = 0. For D^ < D < Dc^, the IH and O^** (where L^^ is replaced by L^^). In the con 



A qualitative phase diagram for (1) is shown in Fig. 1. 

Similar to the spin-1 Haldane phase"^"'"', the spin-2 
IH phase has a bulk gap but exhibits gapless spin-1/2 
edge excitations in open chains. Let us have a closer 
look at how these edge states emerge from Eq. (2). For 
7771 > 0, 7772 < 0, 7713 < 0, our effcctive field theory de- 
scription shares an analogy with Tsvelik's theory-^'' for 
the spin-1 Haldane phase, formulated in terms of three 
Majorana fermions. On a semi-infinite chain, both the- 
ories support three Majorana zero-energy modes at the 
boundary, forming a fractionalized spin-1/2 edge state" '. 
This distinguishes the IH phase from the S0(5) Haldane 
phase (with spin-3/2 edge states^' formed by five Majo- 
rana edge modes ) and the large-D phase (without edge 
states). In fact, we expect that a large family of spin-2 
chains (namely, those described by a similar field theory 
at low energy) also support an IH phase. Interestingly, 
the corresponding Hamiltonians may capture the physics 
of some quasi-lD compounds, and thus the fractional- 
ized edge states in the IH phase may be observed in elec- 
tron spin resonance experiments by doping nonmagnetic 
ions . 

The effective field theory in Eq.(2) also provides the 
order parameters that characterize the three gapped 
phases. To see this, let us view Eq.(2) as five decoupled 
Ising models with the Majorana mass 777 ~ (T — Tc)/Tc, 
where T^ is the Ising critical temperature-'"'" '^. Defin- 
ing the Ising order and disorder operators as Co and 
fia {a — 1,...,5), the S0(5) Haldane, IH, and large- 1? 
phases correspond, respectively, to all five Ising models, 
three Ising models (a = 3, 4, 5), and just one Ising model 
(a = 5) in the ordered phase(s) with (aa) 7^ 0. In order 
to be able to distinguish these phases, we find that it is 
sufficient to use two SOPs 
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FIG. 2: (Color online) String order parameters O^^ and O^". 
They vanish smoothly at D = Dci and D = Dc2 , respectively. 



FIG. 3: (Color online) First 20 coefficients ^^ of the entan- 
glement spectrum for representative points in the (a) SO (5) 
Haldane phase, (b) IH phase, and (c) large-D phase. 



tinuum limit, these SOPs are related to Ising order op- 
erators as ©12 ^ {(Ji){a2),0^^ - (0-3) (c^4) ■ There- 
fore, these SOPs distinguish the S0(5) Haldane (O^^ ^ 
0, ©34 ^ 0), the IH (©12 ^ 0, ©34 ^ 0), and the large-L* 
(0^2 _ (i;^34 _ Q-j p}^a,ses. As we shall see, this is very 
convenient in order to evaluate numerically the phase di- 
agram of the Hamiltonian in Eq. (I). 

The solvability of the S0(5) point D — provides 
an intuitive picture of how these SOPs change with D. 
Starting from Z? == 0, the MPS ground state of Eq. (1) 
has a perfect hidden string order^' : In the S*^ basis, |2) 
and |— 2) appear alternatively in all the configurations 
of the MPS if |0) and |±1) are removed. Similarly, |1) 
and | — 1) also appear alternately, if |0) and |±2) are re- 
moved. This hidden string order is reflected in a nonzero 
value of the SOPs, O^^ = O^"^ = 0.16. When increasing 
D, the uniaxial anisotropy in Eq. (1) tends to suppress 
both the |±2) and |±1) states, but with larger suppres- 
sion strength on |±2). Thus, the string order for the 
I ±2) states is destroyed earlier at 1?^, and an IH phase 
is formed with remaining string order for the |±1) states. 
Here we remind that L^^ and L"^* act on |±2) and |±1), 
respectively. Therefore, we have O^^ = 0, O^^ ^ in the 
IH phase. When increasing D further, the string order 
for the |±1) states is also destroyed at Dca and then the 
system enters the large-D phase with O^^ = O^^ — 0. 

Phase diagram. Let us now explain our numerical 
results for the evaluation of the phase diagram of the 
model. Our technique of choice has been the so-called 
iTEBD algorithm . This algorithm approximates the 
ground-state wave function of the system in the thermo- 
dynamic limit by an MPS. To achieve this, the algorithm 
uses an evolution in imaginary time. The parameter that 
controls the accuracy of the approximation is the size of 
the matrices in the MPS approximation. This parame- 
ter is usually called 'bond dimension', or x- Using this 
method, we have computed MPS approximations to the 
ground state of the Hamiltonian in Eq. (1) for differ- 
ent values of D and x- Then, for each one of these ap- 
proximations we have computed the two SOPs given in 
Eq. (3). In our simulations, we have seen that x = 40 
is enough to reproduce the phase diagram of our model 



with sufficient accuracy for our purposes. 

Our results for the SOPs are shown in Fig. 2. We see 
clearly that the two SOPs distinguish the three phases, 
exactly as predicted by the effective field theory approach 
discussed above. In particular, we obtain the values for 
the critical points Dc, ^ 0.08(1) and Dc^ - 2.63(1). In- 
terestingly, we see that the IH phase actually extends 
through a large region in the phase diagram as compared 
to the S0(5) Haldane phase. 

Entanglement spectrum. The concept of ES was intro- 
duced in Ref.^'' and has proven very useful in the charac- 
terization of ID gapped phases . From a mathematical 
point of view, the ES is simply the spectrum of coeffi- 
cients S,a — — 21ogAQ, where A^ are the Schmidt coeffi- 
cients obtained from the Schmidt decomposition of the 
ground state wave function with respect to a bipartite 
partition in real space, |*) = J2a'^a\^a)\^a)- In t^iis 
equation \'^a ) are the Schmidt vectors for the A{B) 
subsystems. Quite importantly from a numerical per- 
spective, the iTEBD algorithm automatically renders 
this information, since the MPS approximation to the 
ground state wave function is always explicitly written in 
terms of the coefficients Aq, for all possible bipartitions 
of the system into two semi-infinite lines. Thus, the ES 
for these bipartitions can be immediately read out from 
the numerical MPS wave function that approximates the 
ground state. 

In Fig. 3 we show our results for the first 20 coeffi- 
cients <^Q of the ES at several representative points of the 
three phases. As can be seen in the figure, the SO (5) 
Haldane phase is characterized by a quadruple degen- 
eracy in the ES (the coefficients organize themselves in 
quadruplets), and the IH phase by a doubly degenerate 
ES (the coefficients come in duplets). As expected, the 
degeneracy patterns in the ES, introduced by a virtual 
cut, perfectly coincides with the physical edge states in 
these two phases. Also, we see that the large-Z? phase 
has no characteristic degeneracy in the ES. 

The robust degeneracies observed in the ES for the 
S0(5) Haldane and IH phases are protected by the sym- 
metries of the Hamiltonian. For the IH phase, either 



(bond centered) spatial inversion or time reversal sym- 
metry of Eq. (1) is sufficient to protect the doubly de- 
generate ES '. However, the protection of the quadru- 
ply degenerate ES in the S0(5) Haldane phase is be- 
yond the scope of these two symmetries, and is actually 
related to the (Z2 x Z2)^ symmetry. To prove this, it 
is sufficient to show that {Z2 x ^2)^ allows a nontrivial 
four-dimensional irreducible projective representation '. 
Let us identify such a projective representation by fo- 
cusing on the S0(5) point D = 0. To simplify the 
notation, we switch from Cartan basis to vector basis 
for S0(5). Using this notation, the five states in the 
vector basis are written as \n'^) {d = 1, . . . , 5) and the 
S0(5) generators are given by L"'' — i{\n°'){n''\ — h.c). 
The SO (5) point D = has a MPS ground state 
I*) =E{q }Tr(r'^iA---r°"A)|ri''S...,n°"), where A = 
^l4x4 is the (diagonal) matrix of Schmidt coefficients, 
and the four-dimensional matrices F" at each site satisfy 
the Clifford algebra {^^^''} = 2dab^\ Since |*) is in- 
variant under (Z2 x ^2)^ rotations, the local F matrices 
must satisfy the transformation^'^ 



Z^ 



{U" 



)dd' 



^iSab /yab-^^Ydyab 



(4) 



By using the Clifford algebra, we obtain 9ab = 
and the unitaries V'' = ^^[T'',r''], where {t,V^^} x 
{l,y^^} X {1,V^^''} X {1,V^^^} is a four-dimensional ir- 
reducible projective representation of {Z2 x ^2)^, sat- 
isfying {yi^yis} = {^34,^35} = {Fl5,y35| ^ 
and [^1^/34] = [V^^,V^^] = [T/i5,T/34] ^ q. Since 
[V^X] = V y°^ the algebra of V'' guarantees a 
quadruply degenerate ES in the S0(5) Haldane phase. 
Moreover, the (^2X^2)^ symmetry also contains a two- 



dimensional irreducible projective representation, which 
protects the doubly degenerate ES in the IH phase. This 
indicates that the {Z2 x ^2)^ symmetry, besides (bond 
centered) spatial inversion or time reversal symmetry, 
also protects the IH phase. Thus, we conclude that in the 
presence of (Z2 x ^2)^ symmetry, the S0(5) Haldane and 
IH phases are distinct symmetry-protected topological 
phases , which must be separated by topological phase 
transitions (e.g., Dc^ in our model). 

Conclusion. Here we have studied a spin-2 quantum 
chain with an IH phase. By combining effective field 
theory and numerical approaches we have determined the 
phase diagram and the characteristic properties of the 
phases, including edge states and SOPs. An analysis of 
the ES and the {Z2 x ^2)^ symmetry reveals that the IH 
and S0(5) Haldane phases are distinct topological phases 
protected by symmetries. Our results are clear evidence 
of the existence of an IH phase in quantum spin chains. 

Compared to previous work ^ looking for the IH 
phase by considering the spin-2 Heisenberg chain (or 
XXZ chain) with uniaxial anisotropy, we emphasize that 
our starting point of the model (1) is an S0(5) Hal- 
dane phase, which is quite different from the conventional 
Haldane phase in the standard spin-2 Heisenberg chain 
and Affleck-Kennedy-Lieb-Tasaki chain^^'**''. Neverthe- 
less, the IH phase in our model satisfies all characteristic 
features suggested by Oshikawa, and our approach has 
the advantage that the guidance of a low-energy effec- 
tive field theory and full characterization of phases under 
symmetries provide firm evidence of the existence of an 
IH phase. 
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